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We study charmless two-body baryonic B decays in a diagramatic approach. Relations on decay 
amplitudes are obtained. In general there are more than one tree and more than one penguin 
amplitudes. The number of independent amplitudes can be reduced in the large mj limit. It leads 
to more predictive results. Some prominent modes for experimental searches are pointed out. 
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I. INTRODUCTION 

Baryonic modes in B decays are emerging. Decay 
modes, such as B — > D* pn [1], D* pp [2], pp 
ppir [3, 4], Apn~ [5] and A+p [6], have been observed. 
The B° —* A+p decay having 

B(B° -> A+p) = (2.19±°;|| ± 0.32 ± 0.57) x 10~ 5 , (1) 

is the only two-body mode observed so far [6, 7]. By a 
simple scaling of \V u b/V c b\ 2 [8] on the B° — > A+p decay 
rate, the rates of the charmless modes are expected to be 
of the order 10 -7 [6]. This estimation is consistent with 
the (90% confident level) experimental upper limits [7], 



B(B a 



>pp) < 1.2 x 10~ 6 , 
AA) < 1.0 x 10~ 6 , 
B(B~ -> Ap) < 2.2 x 10~ 6 , 



B(B° 



(2) 



obtained by using 31.7 million BB events, and is only one 
order of magnitude below. The number of BB samples 
are accumulating rapidly in B factories. The charmless 
baryonic modes could be just around the corner. 

Motivated by these observations, there are many recent 
theoretical studies in the three-body decay modes [9-17]. 
It is pointed out in Ref. [9] that three-body baryonic 
modes could be enhanced over two-body, by reducing en- 
ergy release to the baryons via emitting a fast recoil me- 
son. The decay rates of D* pn, D* pp, ppK and ApK~ 
modes can be understood to some extent [10-16] and 
the spectra, having threshold enhancement behavior, are 
consistent with predictions [9, 11]. The three-body de- 
cays seem to be much involved than the two-body decays. 
However, in some cases their amplitudes can be related to 
some well measured quantities, such as nucleon magnetic 
form factors, under the factorization approximation, and 
gain better control. 

The two-body baryonic decays are in general non- 
factorizable. One has to resort to model calculations. 
There are pole model [12, 18-20], sum rule [21], diquark 
model [22, 23] and flavor symmetry related [24-29] stud- 
ies. Predictions from various models usually differ a lot 
and early calculations usually give too large rates. Some 
technics developed may still be useful. For example, an 
updated pole model prediction [20] on the B° — > A+p 
rate is consistent with data [6]. 



In this work we use a quark diagram (or topological) 
approach to study charmless two-body baryonic B de- 
cays. This approach was developed and applied to the 
study of the two-body mesonic decays [30-35]. It is 
closely related to the SU(3) flavor symmetry [30, 33]. 
Furthermore, it does not rely on any factorization as- 
sumption. It is stressed that these topological amplitudes 
include long-distant and short distant final state interac- 
tion (FSI) effects [31, 32]. For example, it is used in the 
study on FSI in the B ->■ DP system [35]. 

Motivated by the recent B° — ► A+p observation, the 
topological approach has been applied to the charmful 
baryonic case and the FSI effects are studied [29]. We 
further extend the quark diagram approach to the charm- 
less case and obtain some amplitude relations. In general 
there are more than one tree and more than one penguin 
amplitudes. In principle, these amplitudes can be ex- 
tracted from data. However, so far we do not have any 
relevant data yet. 

It is useful to reduce the number of independent topo- 
logical amplitudes. We use asymptotic relations [36] to 
relate various amplitudes. The same technics has been 
used in the study of the three-body case [12, 14, 15] and 
it leads to encouraging results. For example, the exper- 
iment finding of B(ApTr~) > B(H°p~Tr~) [5] can be un- 
derstood [15] and tree-body decay spectra are consistent 
with the QCD counting rule [37] expectations. Due to 
the large energy release, we expect the asymptotic rela- 
tions to work better in the two-body case than in the 
three-body case. For example, the smallness of two-body 
decay rates may due to some l/m B suppression as ex- 
pected from QCD counting rules. 

The order of this paper is as follows. In Sec. II, we 
formulate the quark diagram approach for the study of 
the charmless two-body baryonic B decays. We consider 
decay modes with decuplet anti-decuplet, octet anti- 
decuplet, decuplet anti-octet and octet anti-octet bary- 
onic final states. Relations on amplitudes are obtained. 
In Sec. Ill, we reduce the number of independent topo- 
logical amplitudes by considering the large ms limit. In 
Sec. IV, we discuss the phenomenology of the charmless 
two-body baryonic decays. We suggest some prominent 
modes for experimental searches. In Sec IV we give dis- 
cussion and conclusion, followed by an appendix for the 
derivation of asymptotic relations. 



2 



II. TOPOLOGICAL AMPLITUDES OF 
CHARMLESS B — > BB DECAYS 

In this section we use a quark diagram (or topological) 
approach to decompose the charmless two-body baryonic 
decay amplitudes. It is useful to re-derive some familiar 
results of the mesonic case first. Since quark diagram 
is a representation of flavor SU(3) symmetry, the topo- 
logical approach should be closely related to the SU(3) 
approach [30]. We use the B — > DP decay to illustrate 
this point and to introduce some useful tools before we 
turn to the charmless case. 

We follow Ref. [38] to decompose B — > DP decay am- 
plitudes according to flavor SU(3) symmetry. We recall 

that the fields annihilating B~ , B d s , creating D°' + , Df 
and creating tt, K, r] S transform respectively as 3, 3 and 
8 under SU(3) [38, 39], 



B = (B~ B° B° s 



D D~ Dz 



( nl + m. 



n 



V 



7T 

K~ 



7T+ 


K+ \ 


tt° , r/s 


K° 


K° 





(3) 



The (du)(cb) operators in the effective Hamiltonian H^j 
can be expressed as {qiH 1 ^) (c6), where q l = (u, d, s) 
and 




// = | 1 




(4) 



It can be understood by using a flavor flow analysis. We 
take the first term of H e g for illustration. In H-w the 
decays are governed by the b — > cq J qi transition with 
the corresponded coupling. The first term of H e ff in 

Eq. (5) is B m D m Hju{. The B m D m part can be in- 
terpreted as a B m to D m transition with the same light 
anti-quark q m (q m ) 1 flavor, while the part is respon- 
sible for the creation of the meson where the IF-emitted 
qiqi pair ends up with. The above picture clearly corre- 
sponds to the external IF-emission topology. Similarly, 
the identification of the C (E) amplitude to the second 
(third) term of H e ff can be understood in the same way. 

It is straightforward to extend the above approach to 
the charmless case and the well known topological de- 
compositions of the charmless B — > PP decay ampli- 
tudes [32-34] can be reproduced. For the b — > uud and 
b — > qqd processes, the tree (Or) and penguin (Op) op- 
erators respectively have the following flavor quantum 
numbers 



Op 



(bu)(ud) = 
{bqiWd) 



(bqi)(q J q k ), 
H k (bQi){q%), 



(7) 



with Hi 2 = 1 = H 2 , otherwise Hf 



H k = 0. The 

j 

flavor structures of |A5| = 1 tree and penguin operators 
can be obtained by replacing d to s and H\ 2 = 1 = H 2 
to Hi 3 = 1 = H 3 in the above expression. By using 
a similar flavor flow analysis as the B — > DP case, we 
obtain 2 



TB m HfU J k Il 



fc n£ir™ + CB m Hf ir^ 



+EB k H) K W l U\ +AB i HfW l W k 



The effective Hamiltonian, in term of the meson degree 
of freedom, for the B — ► DP decay should have the 
same SU(3) transform property of H w . Consequently, 
we have [38] 



H cS = T B m D m R) U{ + C B m Wi 

+EB l mw m D m 1 



H)D J 



(5) 



with probable FSI effects contained in the coefficients. 
The B — ► DP decay amplitudes can be expressed in 
terms of these coefficients [381 



A D o v - 

A D o^o 



= T+C, 



A 



D + 7. 



T + E, 



A D o 



T2 { - c + E) - 



A D+K~ 



E, 



(6) 



The above expression can also be obtained by using 
the topological approach with the coefficients T, C, E 
interpreted as the (color-allowed) external VF-emission, 
(color-suppressed) internal VF-emission and IF-exchange 
tree amplitudes, respectively [31-35]. 

The one-to-one correspondence of the SU(3) parame- 
ters and the topological amplitudes is not a coincidence. 



+PB m H k urul + -PAB k H k U l m IlY\ (8) 

where the A, P and PA terms correspond to annihila- 
tion, penguin and penguin annihilation amplitudes, re- 
spectively. We can reproduce Table I and II of Ref. [33] 
(up to some trivial overall sign changes from wave func- 
tion definitions) by using the above H e g. 

We are now ready to turn to the baryonic cases. We 
will study various decay modes, including dccuplet anti- 
dccuplct, dccuplet anti-octet, octet anti-decuplet and 
octet anti-octet baryonic final states. We start from the 
easiest case, in the sense of flavor structure, and move on 
with increasing complexity. 



A. B to decuplet anti-decuplet baryonic decays 

It is straightforward to extent the quark diagram ap- 
proach to the B to decuplet anti-decuplet decay. As 



1 We use subscript and superscript according to the field conven- 
tion. For example, we assign a subscript (superscript) to the 
initial (final) state anti-quark q m (g m ). 

2 Note that Zfj fc (= H k ) does not lead to any additional term. 
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(c) 



FIG. 1: Pictorial representation of (a) T (tree), (b) P (pen- 
guin) and (c) E (VK-exchange) amplitudes in B to baryon pair 
decays. 



Hence by using the flavor flow analysis and the corre- 
sponding rule, 

qmqi - ©ifci, q l q 3 q m - © Um , (10) 

we have 

H cS = 6T v ^B m H} k V ikl V ljm + 6P v ^B m H k V kil V Hm 

+PA^B k H k V lmn V nml . (11) 

The pre-factors before T v ^ and P v ^ are assigned for 
latter purpose. The above equation is an extension of 
Eq. (8). 

It is useful to discuss the QCD counting rules for these 
amplitudes. In the large tub limit, we need a hard 
gluon to create a qq pair for the tree or penguin topol- 
ogy. An additional gluon is required to kick the spec- 
tator quark in the B meson such that it becomes ener- 
getic in the final baryon pair. In the large tub limit, 
these amplitudes behave like ~ l/m%. The l/m^ fac- 
tor and the Cabibbo-Kobayashi-Maskawa (CKM) coeffi- 
cients lead to a suppressed charmless two-body baryonic 
rates that may be the underlying reason of the nega- 
tive search result. Other topologies, such as VF-exchange 
(Ejyp) ! weak-annihilation and penguin-annihilation 
{PAjyp) topologies, require an additional qq pair creation 
(see, for example Fig. 1(c)) and are further l/m B sup- 
pressed. We should concentrate on T and P amplitudes. 
The decomposition of B to decuplet anti-decuplet de- 
cay amplitudes are given in Table I. We use and 
for the tree and penguin amplitudes, respectively, in 
|A5| = 0(1) processes. 



shown in Fig. 1, we have tree (T), penguin (P) and W- 
exchange (E) amplitudes for processes governed by Ot,p- 
Since there is no external VF-exchange diagram and the 
internal VF-exchange amplitude is not color suppressed, 
we use the symbol T for the tree amplitude. 

The decomposition of decay amplitudes can be ob- 
tained by using the flavor flow analysis. The first step 
is to match the flavor. For example, the decuplet with 
qtq k qi flavor as shown in Fig. 1(a) is produced by the V ik i 
field, while the decuplet with (fqiq" 1 flavor is created by 
the V jlm field, where V jlm is the familiar decuplet field 
with symmetric flavor indices. To be specified, we use 
(see, for example [39]) 



pin 

2)122 
pllS 
2)133 



= A++, V 112 = -^A+ 
V3 



= -=A°, V 22A = A 



1 

— — S* + 2? 123 — - S*° x> 223 — - £*~ 
V3 V6 ' " 

1 =S*°, V 233 = is*-, V 333 = Q~. 



(9) 



TABLE I: Decomposition of B 
tree and penguin amplitudes. 



W amplitudes in terms of 



Mode 


T>T> 


p — 

T>T> 


Mode 


T' — 


1 -D-D 


B~ -> A+A++ 


2^3 


2\/3 


B~ -> E*+A++ 


2\/3 


2^3 


A°A+ 


2 


4 


E*°A+ 




2^2 


e*°iF+ 


V2 


2\[2 




2 


4 


E*-£*° 





2^2 







2V2 







2 







2^3 


A- A 15 





2V3 







2 


B u -> A+A+ 


2 


2 


B" -> E* + A+ 


2 


2 


A ^ 


2 


4 


E*°A° 




2^ 




1 


2 






2^2 







4 


H*~E*^ 





4 







2 







2^ 


A~A^ 





6 


E*-A^ 





2\/3 


B u s -> A+E*+ 


2 


2 


-> E*+E*+ 


2 


2 


A o^o 


v/2 


2^2 




1 


2 


E*°H^ 


V2 


2^2 




2 


4 


E*~E*^ 





4 







4 


z*~TF 





2^3 







6 


A~E*~ 





2^3 


E*~E*^ 





2 



4 



By using the decay amplitudes shown in Table I we we have 
obtain some amplitude relations. For the AS = case, 



A(B~ -» A+A++) - VSA(B- -» A°A+) + A(B~ -» A" A") = 0, 

-» A+A++) = V3 A(B° -> A+A+) = V3A(B° -> A + £*+), 
A(B~ -» A°A+) = -» A AO) = V2 -» A°S^) 

= V2A(B- -» £*°F+) = 2,4(5° -» £*°Fo) = V2 -» S*°S^), 

2>/3A(B _ -» A" A«) = 2A(B° -> A^A - ) = 2\/3^(B° -» A"lF) 

= 3V2 -» E*"!^) = 3 -» E*"IF) = 3.4(5° -> £*~EF) 

= 6A(B~ -> Z*-W°) =6A(B° -> S*~IF) = 2V3A(B° S -» S—TF), (12) 

and for the \AS\ = 1 case, we have 

A(B~ -» £*+A++) - \/6 -» £*°A+) + \/3A(B~ -» £*-A~o) = 0, 

-» £*+A++) = V3^(B° -» £*+A+) = V3 .4(5° -» £*+F+), 

V2 -» £*°A+) = V2A(B° -» £*°A~o) = 2 ,4(5° -» S*°S^) 

= A(B~ -» S*°F+) = V2^(B° -» S*°Fo) = A(5° -» S*°S^o), 

6,4(5" -» S*"A0) = 2v / 3A(5° -» £*~A^) = 6 A(B° S -» E*"XF) 

= 3^2,4(5" -» S*-E^0) = 3,4(5° -» S*~IF) = 3,4(5° -» S^EF) 

= 2^3,4(5- -» fTS^) = 2\/3,4(5° -» fTS^) = 2 ,4(5° -» fi-JF). (13) 



These relations are consistent with Ref. [28]. Note that 
some results in Ref. [28] are obtained by considering dom- 
inant tree or penguin amplitudes only, while the above 
relations include both contributions. 



B. B to decuplet anti-octet, octet anti-decuplet 
baryonic decays 

We extend the previous case to the B — ► BT>, VB cases. 
Note that the (anti-)decuplet parts are as before. For the 
octet part, we use [39] 



strained by 

e ika K + dalK + e akl B- = = ^ b B? + e lbm B{ + e b * m B l b , 

(16) 

which can be shown easily. Hence for each of the qiq k qi 
and (fcpq™ configuration we only need two independent 
terms. 

We are now ready to obtain the effective Hamiltonian 
for the B — ► BV decays. By replacing T>i k i in Eq. (11) 
by tikaB* and e aU Bi , we have 



B = 



E" 



E+ 
_ si _i_ A. 

"0 



n 



(14) 



The has a flavor structure q 3 q a q b e abk - 1 8 j k q c q a q b [39]. 
To match the q%qkQu q l q : 'q m flavor contents of final state 
octet baryons (as shown in Fig. 1), we use 



qiquqi 
q l q 3 q m 



£ikaB l , €ialB k , e a klB i: 



e ljb B, 



m Ibm k?3 



b > 



°6' 



e bjm B 



(15) 



as the corresponding fields in H c g. In fact, not all terms 
in the above equation are independent. They are con- 



5 off = -VQT lB ^B m Hfe ika B^V 1 ^ 
-V6T 2Bl5 B m Hfe akl B a t V 1 ^ 
—V&Pq^ B m H k e kia B l T> hm , 



(17) 



where some pre-factors are introduced for later purpose. 
As argued previously, we shall concentrate on T and P 
terms. Note that we have two tree and one penguin am- 
plitudes. We use T and P' for the |A5| = 1 case. The 
resulting decompositions are shown in Table II. 

There are relations on decay amplitudes. By using 
Table II, we have 
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TABLE II: 


Decomposition of B 


— > 232? amplitudes 


in terms of tree and penguin amplitudes. 




Mode 


T — 

1BX> 


T — 

2BT> 


p _ 

BP 


Mode 


T' 


rpl 

± 1RT> 


p'- 


B~ -> pA++ 


-\/6 


V6 


>/6 


— > E + A++ 


/If 


-V6 




nA+ 


-\/2 





\/2 




V2 





-V2 


AE*+ 


2/V3 


-1/V3 






1/V3 


-i / /o 

1/V3 





E°E*+ 





-1 


-1 




— 1 


1 


2 


E -^o 








-1 










\/2 







o 


V ^ 




n 
u 


u 


1 




-V2 


V2 




1 — > 1 J T— 1-1- A 1 

B — > E A+ 


/n 

V2 


In 

-V2 


-V2 




-72 





\/2 


c°E*° 


1 





-1 


AE*° 


\/2/3 


-I/a/6 


-\/3/2 


AA° 


1/V3 


1/V3 





E°E*° 


o 


-l/\/2 


-1/V2 


E° A 


— 1 


1 


2 


E~E*~ 








-V2 


E~ A~ 


















-V2 


H~E*~ 








\/2 


B" -> »E*+ 




V2 


V2 


B" s -> E + E*+ 


V2 


-V2 


-V2 




-1 





1 


S°H*o 







-V2 


AE*o 


2/V3 


-1/V3 


-V3 


AS^o 


1/V6 


1/V6 





E°E*° 





-1 


-1 


E°E*° 


-1/V2 


1/V2 


V2 


E"S*- 








-V2 


E"E*- 








V2 


s-?F 








-V6 










V2 



TABLE III: Decomposition of B — > VB amplitudes in terms of tree and penguin amplitudes. 



Mode 




T2DB 




Mode 


j" _ 
1™ 


T' - 


p' - 


B~ -> A u p 


V2 





-V2 


B" E*"p 


1 





-1 


E*°E+ 


-1 





1 


H*°E+ 


-V2 





V2 


E*-A 








V3 


H*-A 








V3 


E*"E° 








-1 


H*^EO 








-1 










V2 










a/6 


A"n 








-\/6 


E*"n 








-V2 


B° -> E*°E° 


1/V2 





-1/V2 




1 





-1 


A°n 


V2 


-V2 


V2 


E*°n 


1 


-1 


1 


E*°A 


-1/V6 


n/273 


-^3/2 


S*°A 


-1/V3 


2/V3 


-V3 


A+p 





-V2 




E*+p 





-V2 




E*"E- 








-V2 


H*-E- 








-V2 










-V2 


fi-3- 








-V6 


B" -> E* () H () 


-1 


1 


-1 


-°s — — 


-V2 




-V2 


A°A 


-2/V3 


1/V3 





E*°A 




1/V6 





A°E° 





-1 


2 


E*°E° 





-1/V2 




A+E+ 





>/2 


-V2 


E*+E+ 







-V2 


E* - 3= 








>/2 










y/2 


A"IF 








a/6 


E*-IF 








V2 



V2A(B- -» nA+) + V3^(B" -» AS^+) - -» = 0, 

A(B~ -» pA++) + V3^(B" -» nA+) + 3V2 -» AS^+) - V6 -> = 0, 

A(B~ -» nA+) = -» nAO) = V2A(B2 -» 

-» A(S°)S^+) = V2A(B° -» A(E )S^) = -» A(E°)S ;; o), 

= V3A(B- -» S-S 5 ^) = \/3A(B° -» S-H 177 ) = -» S~?F), (18) 
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for the AS = case and 



V2A(B- -» - V3A(B~ -» AA+) + -» E°A+) = 0, 

A(B~ -» E+A++) - 2\/3 i4(B~ -» + 3V2 -» AA+) - V3A(B~ -» E"A°) = 0, 

i4(fl- -» E+A++) = V3A(B° -» E+A+) = V3A(B2 -» E+E^+), 
A(B~ -» S°E^+) = V2A(B° -» S°E^) = -» S ^), 
A(B~ -» A(E°)A+) = A(B° -» A(S°)A0) = V2A(B° -» A(E°)E^), 
V3A(B" -» XT A°) = E^A - ) = V3A(B° -» E"E^) 

= V6A(B- -» H-E*°) = V3A(S° -» STE* 3 ) = V3A(B° -> H'S* 3 ), (19) 



for the \AS\ = 1 case. 

In Table II we find that the penguin-dominated 
(| AS | = 1) B — > A2? decays do not receive any penguin 
contribution. This can be easily understood. In the pen- 
guin topology, the s quark of A is from the b — > s decay, so 
the u, d quarks of A are correlated to their pair-creating 
partners, the u, d anti-quarks, in the accompanying anti- 
decuplet baryon. The anti-symmetry property of the A 
wave function and the symmetry property of the anti- 
decuplet wave function are responsible for the vanishing 
of the corresponding penguin amplitudes. 



We now turn to the B — ► T>B case. By replacing T) l ^ m 
in Eq. (11) by e ljb B^ and e bjm B[ for the decays, we have 

H cS = -V&T llM B m HfV ikl e lib B^ 
-V&T 2lM B m HfV ikl e bjm B l b 
+V6P vB B m H k V kil e bim B l b . (20) 

We have two tree and one penguin amplitudes. 

The decomposition of decay amplitudes are shown in 
Table III. From the table, we have 



VTiA(B- -» A°p) + V3 A(B° -» A+p) = A(B~ -» A"n) + VTiA(B° -» A°n), 
A(B~ -» A°p) = —\[2A(B~ -» E*°E+) = 2.4(5° -» E*°E°), 
A(B°^A°n) = -V2 A(S° ^ S*°S0), 
A(B°^A+p) = -A(B°^A+S+), 

-» A"n) = -V2 -» E*"A) = V6 -» E*~E°) 

= -V3 -» S*-S0) = V3 A(B° -» S*-S^) = V3A(B° -» E*~lF ) 

= -V3 A(B° -» E*~S^) = -A(B° -» A"E=), (21) 



for AS = processes, and 



V2A(B~ - 


-> E*°p) 


+ A(B° — > 


E*+p) = i4(B~ -» E*~n) + V2,4(B -» E*°n) 


V2A(B~ - 


-» E*°j5) 


= ->1(B- 


-» S*°E+) = V2yl(B° -» S*°S0), 


V2A(B° - 


-» E*°n) 


= 


_> 2*°S°) 


A(B° - 


- E*+p) 


= -^(b2 - 


-» E*+E+), 


V3A(B~ - 


+ E*"n) 







= -\/2A(B- -» H— A) = V6 -» S*-E ) = V3A(B° -» S*~ E~) 

= -V3 A(B° -» S*-S^) = -V3 A(B° -» E*-^ - ), (22) 



for |AS| = 1 processes. The implication of Table II and C. B to octet anti-octet baryonic decays 

Table III on the phenomenology of the corresponding de- 
cay modes will be discussed later. 

The decomposition of B — > BB amplitudes can be 
achieved similarly. To obtain the corresponding decay 
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TABLE IV: Decomposition of B — > BB amplitudes for AS = transitions in terms of tree and penguin amplitudes. 



ivioae 


- 1 IBB 


J 2£B 






p 

IBB 


P 

IBB 


B — ► rip 


-i 

— i 


U 


U 


U 


r 

—0 


U 


AE+ 


-y/2/3 





l/\/6 





-5/\/6 


-5/\/6 


E~ A 














-5/V6 


-5/V6 


E~S° 














-5/V2 


5/V2 


Zj Zj^ 


U 


U 




U 


5/Vz 


— 5/V2 



















—5 


D ► pp 


U 


1 


U 


— 1 


U 


r 




Tin 


1 
1 


1 
1 


n 
u 


n 

U 




c; 


AA 






1 /fi 
1/b 




—0/0 


ZD/0 




1/ V J 


n 

u 


1 /9-/q 


n 

u 


p. /o. /q 
/ z\/ 


c, /9. /q 


Zj 1\ 


u 


u 


i/zyo 


— 1/ V J 


/ Zy 


/9a Fx 
0/ zy -j 


y-y- 

Zj Zj 




n 

u 


n 

u 


n 

U 





p. 


E°E° 








— 1 /2 

l/Z 





—5/2 


5/2 



















5 







-1 





1 





-5 


nA 




-1/V6 








10/\/6 


-5/V6 


nS« 





1/V2 











5/V2 


AH« 






l/\/6 


-1/76 


-5/\/6 


w/VE 


E°H° 








l/\/2 


-1A/2 


5/V2 





E^eF 














-5 






effective Hamiltonian, we replace T>iki (T) 1 ^" 1 ) in Eq. (11) 
by e ika tf and e akl B* {e l ^ b B b n and e«™#£). We have 

H e s = T 1BB B m H l j k e lka B l e ljb B , b n 
+T 2B - E B m Hfe lka B a l € brm B l b 
+T zm B m Hfe akl W^ b BT 
+T iBB B m Hfe ak iB l e bjm B l b 



the eBeB basis, we will have 



-5R 



B m H k e kia B l e hb B r b n 



IBB 



-5 P2BB BrnH e kia B l e tm B b 



(23) 



The coefficients are assigned for later purpose. 

There are four tree amplitudes and two penguin am- 
plitudes. For the tree amplitudes all four combinations 
as suggested in Eqs. (15, 16) are used. For the penguin 
part, only two of the four combinations are independent. 
The two combinations used in the penguin amplitudes 
can be expressed as 

e kia B?e lib B? = \{B,B}% + ^[B,B]f, 



Y2 a bvra yyl 



\{B,B}T 



1 



[B,V\Z - 8?Tr(BB). 



(24) 



One can readily recognize the anti-commutator and 
the commutator parts as D and _F-terms in the usual 
SU(3) combination when dealing with an operator hav- 
ing q k H k n q m flavor quantum number [39]. The Tr(BB) 
term, as denoted as S in Ref. [14], is needed in a non- 
traceless H k case. If we use the D-F-S basis instead of 



^ = ^1+^2), 



Pf = ^(Pi-P2), 



Ps = -P*. 



(25) 

Since we only have two independent components, we have 
a constraint: Pd — Pf + Ps = 0. This is nothing but 
the Okubo-Zweig-Iizuka (OZI) rule in the BB sector. A 
similar constraint is used in Ref. [14], where it is enforced 
by handed. The OZI rule is satisfied automatically in this 
approach, since we are matching the baryon quark flavors 
(see, Fig. 1) in the beginning. 

We shown in Table IV (V) the decay amplitudes for 
IAS 1 ! = (1) processes. From the tables, we obtain 



V3A(B- — > AE+ ) 

V3A(B~ -> E"A) 

A(B°^pp) 
A(B~ -> E°S+) 
A(B~ -» S~W) 

for the AS = case, and 

V3A(B- -» Ap) : 



A(B~ -» S S+) 
+V2A(B- -» rip), 
A(B~ -» T,~TP) 
+V2A(B- -» H-S°), 
-A(B° s ^pZ+),_ 
-V2A(B° -» S°S0), 
-yl(B - S-S^), (26) 



+A(S- -» S°p), 
A(S~ -» S"S0) = >/3A(B- -» S"A) 

+V2 A(S" -» S _ n), 
A(B°^S+p) = -A(S°^S+S+), 



(27) 
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TABLE V: Decomposition of B —> BB amplitudes for |AS| = 1 transitions in terms of tree and penguin amplitudes. 



Mode 


T - 


T' - 


T' - 


T" - 


P' - 


p' - 


B~ -> Ap 


1/V6 





1/V6 





10/V6 


-5/V6 


E°E+ 


-1 











-5 





H"E0 














-5/V2 





E~ A 














-5/V6 


10/V6 


E-n 

















-5 


E°p 







l/\/2 








-5/\/2 


B° -> An 


1/V6 


-1/V6 


1/V6 


-1/V6 


10/V6 


-5/^6 


S"E- 














-5 





E°E° 













5/v^ 





H°A 


-1/V6 


VW3 








-5/\/6 


10/V6 


E+p 





-1 





1 





-5 


E°n 


-1/V2 


1/V2 


1/V2 


-i/\/2 





5/V2 


AA 


-1/3 


1/6 


-1/3 


1/6 


-10/3 


5/3 




-1 


1 








-5 


5 
















-5 


5 


E+E+ 





1 





-1 





5 


E°E0 





1/2 





-1/2 





5 


E"EF 

















5 


AE° 





-l/2y/3 





-1/2^3 








E°A 


1/V3 


-l/2y/3 


-1/V3 


l/2y/3 









for the IAS 1 ! = 1 case. More relations in the |A5| = 1 contribution. We have 
case can be obtained by neglecting the sub-leading tree 



A(B~ - 

V2A(B- 
A(B~ - 
A(B~ - 



+ H~A) = A(B°->Z°A), 

S-SF) = A(B° S -+ Z°ZP), 

-» Ap) = V2 A(B a -» An) = V3 A(B° S 

S°E+) = V2A(B- -» S-S°) = ^1(B 

■* E-n) = \/~2A(B~ -» E°p) = - 
= -» E+E+) = -» 



-AA), 

-» S-S - ) = -V2 -» S°E°), 
» E+p) = -V2 -» E°n) 

E°E°) = -» E~E~ ). 



(28) 



Rate relations implied by Eq. (28) are consistent with 
those obtained in Rcf. [27]. In fact, Rcf. [27] uses a 
generic SU(3) analysis and has three independent pen- 
guin amplitudes. One can reduce these amplitudes into 
two penguin amplitudes by imposing the OZI rule. As 
noted before the quark diagram approach is consistent 
with the OZI rule (see Eq. (25)). As a result, we have 
more relations. 

To summarize, based solely on the flavor structure of 
i/wj we obtain one tree and one penguin amplitudes for 
the DT> case; two tree and one penguin amplitudes for 
each of the BT> and the BT> case; four tree and two pen- 
guin amplitudes for the BB case. In principle, these am- 
plitudes can be extracted from data. Various relations 
on decay rates as implied in Eqs. (12, 13, 18, 19, 21, 
22, 26, 27, 28) are obtained and can be checked exper- 



imentally. It should be pointed out that although we 
concentrate in tree and penguin amplitudes, the exten- 
sion to include other contributions, such as M^-exchange, 
electroweak penguin, weak annihilation and penguin an- 
nihilation amplitudes, is straightforward. 



III. REDUCTION OF TOPOLOGICAL 
AMPLITUDES 

Up to now we only make use of the flavor structure 
of the weak Hamiltonian. Amplitudes are decomposed 
topologically. In most cases, we need more than one 
tree and more than one penguin amplitudes. Although 
some relations are obtained, it is a long term experimen- 
tal project to verify them. In fact, by taking into account 
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A(B - 


- B 1 B 2 ) 


A(B- 


-» 2?iB 2 ) 


A(B- 


-» Bi© 2 ) 


A(B - 





of the chirality structure of i?w 5 the number of indepen- 
dent amplitudes is reduced in the large mj limit. It leads 
to more predictive results. 

In general the charmless two-body decay amplitudes 
can be expressed as [12, 19] 



= u i( A bb + ^5B b1S )v 2 , 

= i <(A DB + 7 5 .& D ;g)w 2 , 

TUB 
qV- 

H —u^Cjyp + j 5 D V v)v2vl29) 



where q = p\ — P2 and u^, are the Rarita-Schwingcr 
vector spinors for a spin-| particle. The vector spinors 
in various helicity states can be expressed as [40] 
«„(±§) = e„(±l)«(±|), u„(±|) = ( c „(±1)u(t§) + 
\pl e M (0)u(±5))/v / 3, where e M (A) and u(s) are the usual 
polarization vector and spinor, respectively. By us- 
ing q ■ e(A)i,2 = T&\fimBPc/mi,2, where p c is the 
baryon momentum in the B rest frame, and the fact that 
e i(0) • e 2 (0) = (r/ig — m\ — to|)/2toiTO 2 is the dominant 
term among e^(Ai) • e 2 (A 2 ), we have 



decay as an example. According to the above argu- 
ment, the final state (anti-)baryons should have u^d^q^ 
{q' L URSL) configuration from the tree b — ► ulurcIl decay 
and d L q L(R) q' R(L) (q' L ( R )q R ( L )S L ) configuration from the 
penguin b — > dLqL( R )Q R (L) decay. In the large tub limit, 
as the spinor helicity identify to chirality, we should have 
B — > —A in the above equation. 

The number of independent topological amplitudes is 
reduced in the large tub limit. The corresponding asymp- 
totic relations are derived in the appendix. Some of the 
reduction can be understood by using the chirality struc- 
ture of weak interaction, the helicity argument and the 
anti-symmetry property of the A wave function. 

The b — > ULUftdL process have final state u and d 
quarks both with left-handed chirality or helicity (as light 
quark masses are neglected). It is well known that in the 
A wave function the u and d quarks have opposite helicity 
(as shown in Eq. (A2)). Therefore, the tree amplitudes 
in the AS = 0, B — > AV, AB decays should be vanishing. 
As shown in Table II and Table IV these amplitudes are 
proportional to 2 T 1B ^ — T 2B ^ and 2T 1 
T iBB , respectively. We should have 



3BB' 2T 2SS~ 



T (')_ _ 
1BV 



i T c)_ 

2 2BV 1 



(32) 



and 



A(B -> V X B 2 ) 
A(B -» B 1 V 2 ) 
A(B -» Di© 2 ) 



2 p c 



t o — ( A vb + 75S P g)w 2 , 

O TOl 

i \f\^ 2 U ^ A BV + 75-B B5 )u 2 , 

9 

^-, l( ^ + 75 ^), 2 , (30) 



where A^— = - 2{p c /m B fC v ^ 1 B'^ = B v ^ - 

2(p c /rriB) 2 D T yp and decuplets are in h = ±| helicity 
states. Note that the first two equations are exact while 

3> m\ 2 , TOiTO 2 limit. For 
'2 _ /■„ , ^\ 2 



the last one holds in the m 2 B 



these amplitudes to have the same m B = (pi + P2) be- 
havior in the large m 2 B limit, the extra factor of p c and 
itlb in Eq. (30) should be suitably compensated by the 
l/m B power of the corresponding A and B terms. For 
all four B -» BiB 2 (BB = BB, VB, BV, VV) decays, 
we can effectively use 



A(B -» BiB 2 ) - Si (A + 75S)v 2 , 



(31) 



T (')_ = i_ T w_ 

IBB 2 3Ke ' 



T (')_ = i_ T (')_ 
2BB 2 ABB ' 



(33) 



Although above relations are obtained in the AS = 
shown in the appendix, they also hold in the 
\AS\ = 1 case. 

Furthermore, we note that there should be no penguin 
contribution to the B° — > AA decay amplitude. In the 
penguin type b — > dt decay, the quark content of the 
produced A should be diiijjSi, since the u d pair of the 
A wave function is in a spin zero configuration. From the 
helicity (or chirality) conservation, the correlated anti- 
baryon is in a left-handed helicity state with quark con- 
tent dLULSR, which cannot match to the wave function 
of A and resulting a vanishing penguin contribution on 
the B° — ► AA mode. As shown in Table IV, its penguin 
amplitude is proportional to P 1BB 
have 



5P 2BB and we should 



P 



IBB 



5R 



(') 



2BB' 



(34) 



as the decay amplitudes. 

The chiral structure of weak interaction provide further 
information on A and B. For example, in the AS = 
processes, we have b — ► ULU R dL and b — > rfL<7L(fl)te(L) 
decays. The produced c?l quark is left-handed in both 
tree and penguin decays. Furthermore, as strong inter- 
action is chirality conserving, the pop up quark pair q'q' 
should be q' L q' R or q' R q' L - From the conservation of helic- 
ity, the produced baryon and anti-baryon must be in left- 
helicity states. To be more specify, we take the B° s — > BB 



According to the decomposition shown in Eq. (25), this 
relation can be expressed by (Pd/Pf) —> (3/2), which 
is consistent with the asymptotic relation on the baryon 
scalar and pseudoscalar form factors obtained in Ref. [15]. 

There is no constraint from processes involving A final 
states in the B — > VB decays and we still have two inde- 
pendent tree amplitudes in B — > BB decays. Asymptotic 
relations lead to further reduction, 



r('_ 

1VB 



r (')_ 
2T>B' 



(35) 
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and 



p(')_ 

IBB 



r (')_ 
2BB' 



(36) 



in the large mg limit. 

It is interesting to note that in Ref. [36] the octet-octet 
and octet-decuplet systems are related asymptotically. 
Similarly, we have 



T (>) = 
p(') 



r (0_ 

WD 



J>(')_ _ J>(')_ _ Ji(' 



p(')_ 

vv 



' IBV 

P (0_ 



"1£>B 
p(')_ 



IBB 

p(')_ 
IBB' 



(37) 



in the large wis limit. In that limit, we need only one 
tree and one penguin amplitudes under the quark dia- 
gram approach for all four classes of charmless two-body 
baryonic modes. In the next section, we will use the 
above equation as an approximation to study the decay 
rates. It should be verified by data that whether this is 
a good approximation or not. 



IV. PHENOMENOLOGICAL DISCUSSION 

In this section, we discuss the phenomenology of 
charmless two-body baryonic modes by using the formal- 
ism developed in previous sections. Since none of the 
charmless two-body baryonic decay is observed so far, we 
suggest some prominent modes for experimental searches. 
This section is divided into two parts, discussing tree- 
dominated and penguin-dominated modes, respectively. 
We summarize our results at the end of this section. 



A. Tree-dominated modes 

For AS — modes, we expect tree amplitudes to dom- 
inate. We can estimate their relative rates by neglect- 
ing penguin contribution. Rates are normalized to the 
B° — > pp rate. As noted before, a simple scaling of 
\V u b/V c b\ 2 on the B° — ► A+p decay rate hints at a <~ 10~ 7 
rate for the charmless case [6]. A pole model calculation 
also gives B(B° — ► pp) = 1.1 x 10~ 7 [12]. For illustration, 
we take B(B° — > pp) = 1 x 10~ 7 as the reference rate for 
these tree-dominated decay rates. 

As noted in the end of the previous section, we use 
the asymptotic relations stated in Eqs. (32-37) for tree 
and penguin amplitudes shown in Tables I-IV. Hadron 
masses and B~'°, B s lifetimes taken from Ref. [8] are used 
to indicate some SU(3) breaking effects. Results obtained 
are given in Table VI. For comparison, we quote results 
of pole model [12], diquark mode [23] and sum rule [21] 
calculations in the subsiding order in parentheses. We 
normalize the tree contributed Bt(B° — > pp) of Ref. [12] 
and Ref. [21] to 1 x 10" 7 . Note that the B° -» pp rate 
only decreased by 16% by including penguin contribu- 
tions in Ref. [23]. The tree-penguin interference effect is 
not prominent. We will wait for data before carry out a 
detail study of the tree-penguin interference effect. 



Rate relations shown in Table VI are consistent with 
Eqs. (12J8, 21, 22, 26). Wc find that the tree-dominated 
B — > NA modes have the structure 



B(B- 



P A++) ~ 3B(B° 
~ 3B(B~ 



P A+) 

>nA+)~3S(B"^nA0), 
(38) 



which is consistent with Refs. [12, 19], while the B(B — ► 
NA)/B(B° — > pp) ratios are smaller than those obtained 
in Ref. [12] by a factor of two. We have B(B° — > pp) ~ 
B(B~ — > pn), which is close to the sum-rule result [21]. 
On the other hand, our prediction on B(B° — > nn) : 
B(B° — > pp) ~ 4 : 1 is different from all quoted earlier 
results. This different may be related to the following 
fact. It is easy to see from Table IV and Eqs. (33, 36) 
that by neglecting the sub-leading penguin contribution 
we have 

A(B~ -» np) - A(B° -» pp) = A(B° -» nn). (39) 

The origin of the above equation is the correspond- 
ing er coefficients shown in Table IX. Note that the 
above equation is different from the AI = ^ rule re- 
lation [12, 19], which implies \A(B° -> pp) - A(B a -> 
nn)\ = \A(B~ — > np)\. In fact, a highly suppressed 
B(B~ -» TT-Tr a )/B(B° -> 7T+7T-) ratio is followed from 
the AI = i rule in the charmless mesonic sector. Al- 
though the present data is still fluctuating, it suggests 
B(ir~Tr°)/B(Tr + Tr~) ~ 1 [41]. We do not worry about the 
deviation of the AI = | rule in Eq. (39) at this stage. 
Note that we have B(B° — > AA) = 0. As shown in the 
previous section, the vanishing of the B° — > AA rate is 
due to the anti-symmetric property of the A wave func- 
tion and the chirality argument in the large mg limit. 
This mode provides a useful check of the approximation 
made. 

Upper limits on B° — > pp, AA and B~ — > Ap are 
obtained [7] (see, Eq. (2)). Some observed three-body 
modes [3-5] also provide constraints. To discuss the fea- 
sibility of experimental observation, it is useful to recall 
that B{YP -» A 7 ) ~ 1, B(E -» Att) - 1, B(A+ -» pn ) ~ 
§, B(A° -» pn') ~ i and £(£* -» Att) = 88% [8]. 

As shown in the first column of Table VI, there are ten 
B° decay modes having rates larger than the reference 
B° — > pp rate. We can search for the E°A mode, which 
has a rate about three times of the pp rate and clear sig- 
nature for reconstruction. Up to now, most experimental 
searches in baryonic modes do not involve n° in recon- 
struction. With o ne n ° in the final states, we can search 
forpA+, A + p, and S*°A modes, which have rates 

about twice of the B° — > pp rate. Note that B a — > nn, 
A°n decays have rates larger than the pp rate by factors 
of four and eight, respectively, but require n or even n in 
detection. 

We now turn to decay modes as shown in the 
second column of Table VI. One should search for the 
A+A++ mode, which have the largest rate, B(B~ — > 
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TABLE VI: Decay rates for AS = tree-dominated modes. We consider only tree amplitude contribution. Rates are normalized 
to B(B° -> pp). For comparison, we show results of pole model [12], diquark mode [23] and sum rule [21] calculations in the 
subsiding order in parentheses. 



Mode 




B{10- ? ) 


Mode 




B(W~ 7 ) 


Mode 


Z3(10- 7 ) 


B° -> pp 


1 a 


(1.1, 0.84, 1) 


B~ — > np 


1.09 


(5.0, 0, 0.6) 


B l > -> pE+ 


0.96 


nn 


4.00 


(1.2, 0.84, 0.3) 


E°E+ 


1.99 




nA 


1.46 


AA 





(0, 0.38, — ) 


P A++ 


6.19 


(14, 0.63, 0.25) 


nE° 


0.48 


E°A 


2.79 




nA+ 


2.06 


(4.6, 0.70, — ) 


E°H° 


7.17 


E°E° 


0.91 




T,°W+ 


3.81 




pW+ 


1.85 


pA+ 


1.90 


(4.3, 0.28, 0.1) 


A°p 


2.06 




nE*° 


0.92 


nA° 


1.90 


(4.3,0.28, -) 


E*°E+ 


0.95 






3.40 


E°E*° 


1.75 




A+A++ 


11.72 




A+E+ 


1.83 


A+p 


1.90 




A°A+ 


3.91 




A°A 


2.78 


A°n 


7.60 






1.82 




A°E° 


0.91 


E*°A 


1.34 












3.44 


E*°E° 


0.44 










a + tF+ 


3.50 


A+A+ 


3.60 










A°E*o 


1.75 


A (, A^ 


3.60 












1.63 



E*°E*° 0.84 



°We take B(B° — > pp) = 1 X 10 7 as a reference rate. 



A+A++) = 11.7 B(B° -> pp), and only reduces ~ 30% 
in producing p7r°jj7r~ final state. Two modes, pA ++ 
and A°p, can decay to the all charged ppir~ final state. 
They have rates six and two times of B(B a — > pp) and are 
searchable. Their rates are within the tree-body (total) 
rate constraint, B(ppn-) = (3.06±^ ± 0.40) x 10~ 6 [4], 
and should be accessible. 

For the B s case, one needs all charge track for de- 
tection. As shown in the third column of Table VI the 
and A°A modes have rates of one to three times 
of B(pp) and should be searchable. 

There are some pure-tree modes in IAS 1 ! = 1 pro- 
cesses. We use T' — TV* s /V* d to estimate their rates 
and obtain B{B~ -» AA+) = 0.15 x 10~ 7 and B(B° S -» 
AS*°, S*°A) = 0.07 x lO" 7 , which arc quite small as 
expected. There are some pure-penguin modes in the 
AS = process as shown in Tables I-IV. Their rates 
are suppressed and will be briefly discussed in the next 
sub-section. 



B. Penguin-dominated modes 

The relative rates for penguin-dominated modes can be 
obtained similarly as in the previous subsection, but in- 
stead of neglecting penguin amplitudes we neglect tree 
amplitudes. We use B~ — > Ap as a reference mode. 
The present experimental limit is B(B~ — > Ap) < 
2.2 x 10~ 6 [7]. On the theoretical side, a recent pole 
model calculation gives B(B~ — > Ap) = 2B(B° — » 
pp) = 2.2 x 10 -7 [12], while a diquark model calcula- 
tion gives B(B~ -» Ap) = 0.18 B(B° -» pp) [23]. Due 
to the controversial situation of these predictions, we use 
B(B~ — > Ap) = lxl0~ 7 asa convenient reference rate 



for illustration. In Table VII we show rates for penguin- 
dominated modes, while rates for pure-penguin modes 
are given in Table VIII. Rate ratios are consistent with 
Eqs. (13, 19, 22, 27, 28). For comparison, we show results 
of pole model [12], diquark mode [23] and sum rule [21] 
calculations, in the subsiding order in parentheses. The 
diquark and sum rule results are normalized as before. 

The B(B° An)/B(B- -> Ap) ratio is similar to 
those obtained in pole model and diquark model calcu- 
lations [12, 23], while the B(W -> Y,+p)/B(B~ -> Ap) 
ratio is consistent with the pole model result [12]. For 
other cases, we have quite different predictions. 

As shown in Table VII, the B~ — > Ap decay is still 
the best mode to search for. Other modes may have 
much smaller rates or much lower detection efficiencies. 
For example, although the B~ — ► S°p decay signature is 
easy to identify, its rate is too small to search for. This 
smallness as compared to t he A p rate can be traced to 
the e P ratio (-1/3V2 : ^^3/2) as shown in Table IX. 
Note that the same ratio is obtained in the Ap vs. E°p 
(pseudo) scalar form factors in the asymptotic limit [15]. 
The resulting prediction , B(Ap~Tr~) > B(Y,°pir~), is sup- 
ported by data [5]. The B~ — ► rate is about 
twice the Ap rate, but we need two ir° in reconstruc- 
tion. Note that an all charged final state can be found 
in the B~ -> £*+A++ -» Att+Jjtt- decay with only 10% 
reduction in rate. 

The penguin-dominated B° decay rates are shown in 
the second column of Table VII. In general t heir rate 
are within B(B~ -> Ap). We note that the S°A" final 
state has one third of chance to decay to YPpir~ and 
the rate is within the three-body total rate constraint 
{B(Y?piT-) < 3.8 x 10^ at 90% confidence level [5]). 
For the B s case, the AA mode is the most searchable 
one. 
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TABLE VII: Decay rates for |AS| = 1 penguin-dominated modes. We consider only penguin amplitude contribution. Rates 
are normalize to B(B~ — » Ap). For comparison, we show results of pole model [12], diquark mode [23] and sum rule [21] 
calculations in the subsiding order in parentheses. 



Mode 




B(W- 7 ) 


Mode 






Mode 


Z3(10- 7 ) 


B~ -> Ap 


1 b 


(2.2, 0.18, < 3.8) 


W E+p 


0.07 


(0.2, 0.88, 7.5) 


B u s -> E+E+ 


0.06 


E°p 


0.04 


(0.6, 0.18, 3.8) 


An 


0.92 c 


(2.1, 0.21, — ) 


AA 


0.60 


E°E+ 


1.70 




E°n 


0.03 c 


( — , 0.21, — ) 


E°E° 


0.06 


E°A+ 


0.28 


( — ,0.08, — ) 


H°E^ 


0.78 






0.98 


H°E*+ 


0.13 




E o A0 


0.26 


( — 17 — ) 


E+E*+ 


0.12 


E+A++ 


0.42 


(2.0, 1.10, 7.5) 




0.06 




E°E*° 


0.12 


E*°p 


0.07 




E+A+ 


0.13 


(0.6, 0.57, 7.5) 




0.12 


H*°E+ 


0.13 




E*°n 


0.06 




E* E« 


0.12 


E*°A+ 


0.53 




E*+p 


0.13 




H*°H° 


0.12 


E*+A++ 


0.80 




H*°A 


0.20 




E*+E+ 


0.12 


H*°E*+ 


0.98 




2*°eo" 


0.06 






0.88 








E*°A° 


0.49 




E*°I>° 


0.24 








E*+A+ 


0.24 




E*+E*+ 


0.24 



H*°E*° 0.45 



''We take B(B~ -<• Ap) = 1 x 1CT 7 as our reference rate. 
These modes may have large tree contributions. See text for 
detail discussion. 



TABLE VIII: Rates (normalized to B(B~ -> Ap) = 1 x 10~ 7 ) for pure penguin modes in the |AS| = 1 process. For comparison 
we show results of Ref. [12] in the parenthesis. 



Mode 


B(10- 7 ) 


Mode 


B(10- 7 ) 


Mode 


B{10- 7 ) 


B~ -> H~E° 


0.85 


B° -> H-E^ 


1.56 




0.98 


S"A 


0.10 


H°A 


0.10 


E"E- 


0.06 


E"n 


0.07 


S"E*- 


0.12 


E-E^ 


0.12 


H~E*° 


0.06 


E~A~ 


0.38 




0.12 


E^AO 


0.14 (0.7) 


E*-lF 


0.12 


E'-fF 


0.12 


E*-n 


0.14 


Q~E- 


0.33 




0.12 


H*"A 


0.20 




0.73 


E*~E*^ 


0.24 


H*-E° 


0.06 


E*~E*^ 


0.90 




0.88 


n~E° 


0.36 




0.62 




1.82 


E* _ A° 


0.27 










^"*~ E*° 


0.49 












0.68 











We note that some modes may have large tree contri- 
butions. For example, as one can see from Table V and 
Eqs. (33, 36) that B° — > An, T,°n amplitudes have rel- 
atively large T coefficients. By using V = TV*JV* d , 
we have B T (B° -> An, £°n) = (0.12, 0.07) B(B° -> pp) 
from tree contributions. For these modes, tree-penguin 
interference may affects their rates dramatically. For ex- 
ample, by including tree contribution in B° — > E°n de- 
cay, Ref. [23] obtains a factor of five enhancement of an 
earlier calculation [22]. As one can see from Tables I-V 
and asymptotic relations that other IAS*! = 1 modes do 
not have such a large relative tree contribution. The tree- 
penguin interference for them are expected to be much 
milder. 

Rates for pure-penguin modes are shown in Table VIII. 
They are in general not favorable than the B~ — > Ap 



mode. The largest mode of pure-penguin B~ decays is 
S~£° and one may search for it after the observation of 
the Ap mode. Although the B° — > S~£~ rate is 1.6 times 
of Ap rate, we need A7r _ 7T7r + for detection. One may 
have chance to search for B® — > S~S~, decays 
through A7r - A7r + , AK~AK + final states, respectively. 

As noted in the end of the previous subsection, there 
are some pure penguin decays in AS* = processes. We 
can estimate their rate by using P = P' V t * d /V t * s . Their 
rates are ty pica lly of 10~ 9 and even the largest rate, 
B(B° -> A L A-) ~ 9 x 10~ 9 , is still undetectable. 

We give a summary of our suggestions before ending 
this section. We find that in AS = processes, in ad- 
dition of the B° — > pp search, it is useful to search for 
the B° -> S°A decay and B~ -> A +A++, pA++, A°p 
decays. In particular, the A + A++ rate is larger than 
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the pp rate by one order of magnitude. It could be the 
first charmless two-body baryonic mode to appear. For 
\AS\ — 1 processes, the B~ — > Ap decay is still the 
best mode to search for. After the observation of any 
of the above mention modes, one should also search for 
other sub-dominated modes, such as B° — > pA+, A + p, 
s *o-£ ; s _jp and B _ s *+ a+T decays. 

For the B s case, as all charged tracks final states are re- 
quired for detection, one can search for the p£(*)+, A°A, 
and modes, which have rates compatible 

to B° — > or — > Ap decay rates. We note that 
although rates shown in Table VI seems more promising 
than those shown Table VII and Table VIII, the refer- 
ence B~ — ► Ap rate may be greater than the reference 
B° — > pp rate and resulting a different prospect. Never- 
theless, there are many promising modes, as good as the 
pp mode, in AS = processes, while the Ap mode is still 
be best candidate in |A5| = 1 processes. 



V. DISCUSSION AND CONCLUSION 

In this work, we use a quark diagram approach to clas- 
sify various contributions on charmless two-body bary- 
onic B decays. This approach is closely related to flavor 
symmetry, since a quark diagram is a representation of 
it. Relations on decay amplitudes are obtained. As noted 
in Sec. II, our results are consistent with some generic 
SU(3) results with the OZI rule imposed. Since we match 
baryon and qqq flavors in the beginning, the OZI rule 
should be satisfied. In this sense, we may view the dia- 
gramatic approach as an flavor symmetry approach with 
the OZI-rule constraint. 

We concentrated on tree and strong penguin contribu- 
tions in decay rates. Other contributions are sub-leading 
and can be included later. For example, the electroweak 
penguin contribution can be included by a simple re- 
definition of tree and strong penguin amplitudes as in 
the mesonic case [34]. Some contributions, such as weak 
exchange and annihilation topologies, can be included 
by extending Eq. (11) to other cases, but they are ex- 
pected to be l/m B suppressed. Amplitude relations are 
obtained by considering tree and penguin contributions. 

There are in general more than one tree and more than 
one penguin amplitudes to describe charmless two-body 
baryonic decays. Reduction in the number of indepen- 
dent amplitudes can be achieved by considering the chi- 
rality nature of weak interaction and the large tub limit. 
Some relations can be understood by using the anti- 
symmetry behavior of the A wave function. An interest- 
ing check is the experimental verification of the vanish- 
ing or highly suppressed B° — > AA rate. As in Ref. [36], 
the octet-octet and octet-decuplet systems are related 
asymptotically, we need only one tree and one penguin 
amplitudes to describe charmless two-body baryonic de- 
cays. The followed results are very predictive and can be 
checked by experiments. 

We suggest some prominent modes to search for. In 



AS = processes, there are many promising modes. In 
addition of the B° — > pp search, it is useful to search 
for B° -> S°A decay and B~ -> A+A++, pA++, A°p 
decays. For \AS\ = 1 processes, the B~ — ► Ap decay 
is still the best mode to search for. After the obser- 
vation of any of the above mention modes, one should 
also search for other sub-dominated modes. For the B s 
case, one can search for the p£(*)+, A A, S~S~ and 
ri 17 decay modes. Although the above suggestions are 
obtained by considering the dominant contributions, we 
do not expect large modification of relative rate ratio in 
most cases. The tree-penguin interference effects can be 
included later after the appearance of data. 

To conclude, we use a quark diagram approach to study 
charmless two-body baryonic decays. The topological 
amplitudes can be extracted from data. We further ap- 
ply asymptotic relations to reduce the number of inde- 
pendent topological amplitudes and obtained predictive 
results. We have pointed out several promising modes 
to be added to the present experimental searching list. 
The discovery of any one of them should be followed by 
a bunch of other modes. 
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APPENDIX A: ASYMPTOTIC RELATIONS IN 
THE LARGE m B LIMIT 

We follow Ref. [15, 36] to obtain the asymptotic rela- 
tions stated in Sec. III. As noted we only need to consider 
helicity ±i states. The wave function of a right-handed 
(helicity= i) baryon can be expressed as 

|B ; T) ~ -±={\B ; TIT) + |B ; HI) + |B ; ||T)), (Al) 

i.e. composed of 13-, 12- and 23-symmetric terms, re- 
spectively. For B = A, p, n, S°, A, we have 

|A++;TIT) = u(l)«(2)«(3)| TIT), 
|A-;TIT) = d(l)d(2)d(3)| TIT), 

|A+;TIT) = -^[«(l)«(2)d(3)+«(l)d(2)u(3) 

+d(l)u(2)u(3)]| TIT), 
|A°;TIT) = (|A+;TIT) with u <-» d), 
|X*+;TIT) = (|A+;TIT) withers), 

|£*°;UT) = -7=Ml)d(2)a(3)+ permutation] | TIT), 
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Ip;UT> 

|n;UT> 
|£°;UT) 



u(2) 



d(l)u(3) +u(l)d(3) 

-y|«(i)d(2)«(3) | TIT), 
(— |p ; TiT) with u d), 

u(l)d(3) +d(l)u(3) 
u(2)d(3) + d(2)u(3) 
u(l)d(2) + d(l)u(2) 



*(2) 



+ 



5(1) 



+ 



|A;TTT> = 



2V3 

d(2)u(3) - u(2)d(3) 



a(3) 
S (l) 



| u(l)d(2) - d(l)n(2) 



TIT), 



TIT),(A2) 



for the corresponding |B ; TIT) parts, while the 12- and 
23-symmetric parts can be obtained by permutation. To 
be consistent with Eq. (14), our A state has an overall 
negative sign with respect to that of Ref. [36] . 

Following Ref. [36] and using the helicity argument in 
Sec. IV, asymptotically we have 



q R 



q R 



q R 



*L(R) 



lR(L) 



q L b 



(a) 



q L 



i 



(b) 



d L , S] 



q R 



d L> S L 
q L(R) 

q R (D 



(B(p 1 )|0|BB , (p 2 ))=u( P i) 



1-75 



F(t) 



F(t)= E ei(B'-B-B)Fi(t), (A3) 

i=T,P L ,P„ 



FIG. 2: (a) Tree and (b) penguin B' — B — B diagrams in the 
asymptotic limit. 



where O are the operators in H eS . For simplicity, we TABLE IX: The coefficients e T ,p(B' - B - B) for various 
illustrate with the space-like case. Note that the above modes obtained from Eqs. (A4, A5). 
equation is obtained in the large t(= (pi — P2) 2 ) limit, 
where we may take a large tub limit or view the B meson 
as a virtual particle with a large momentum. Quark mass 
dependent terms behave like m q / -y/\t\ and are neglected. 
As shown in Fig. 2(a), the tree operator (ub)v-A(du)v-A 
can generate a ~B'(q' R Ul qR)~B(q' L b)-B(uL 4l qn) cou- 
pling. The corresponding coefficient ex(B' — B — B) is 
given by 



e r (B' -B-B) 

= (B; ||T \Q[q' R (l) - «i(l); u L (2) - d L (2)]|B' ; TIT) 



B' - B - 


B 


ex 


ep 


B' 


~ B - 


B 


eT 


ep 


E* u - B u - 


-E*° 


1 

3 


2 
3 


A°- 


B~ - 


- A" 





2 

V3_ 
2V2 
3 


A+ -B- - 


-A 


2 
3 


4 
3 


A+- 


- S" - 


-E*° 


V2 
3 


A++-B- 
A+-B" 


-P 
-P 


V2 \/2 


E*° 
A+ 


-B° 
-B~ 


- A 

— n 




v/2 
3 


1 

n/6 
\/2 
3 


E*+-B l > 


-P 






A+- 


- B~ 


-E° 


1 


2 


p-B~ - 


A 


3 




P- 


B°- 


A+ 


3 


3 


n-B~ - 
p-B~ - 


A" 

E*° 




1 

-!- 


-7i 

1 

-T 2 - 


A- 

n — 


B°- 
B°- 


E*° 
E*° 


1 

s/6 
2 


1 
s/6 
1 

"A" 



+<B;TU|Q[9k(3)^«L(3);« L (2)^d i (2)]|B';UT>, "-if-n -§ -§ A-B°-E° ^ 

(A4) S°-B°-E° -i -I E°-B°-A j. 

A-B°-A n-B°-A ^1 \/l 
where Q[^(l(3)) -» Ui (l, 3); u L (2) -» d L (2)] changes p-B~-TP p-B~-A 
the g'(l(3))| T>®«(2)| |) part of |B'; TIT) to the u(l(3))| | = 
) ®d(2)| |) part. 

Similarly coefficients ep Lt p R (B' — B — B) for 
theB g L <&)-B(g' L b)-B(d L q L q' R ) and B'(^ q R q'[)- 
B(q' L b)-B(dLqRq'l) couplings governed respectively by 
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the penguin operators {db)v-A{qq)vTA & re given by 
e PL (B'-B-B) 

= (B; lit \QWr(1) - d L (l); q L (2) -> <fc(2)]|B' ; TIT) 

+ (B; TU \QW R (3) - d L (3)^ L (2) -> gi (2)]|B' ; TIT), 
e P (B' - B - B) = e PL (B' - B - B) 

= e PR (B' -B-B). (A5) 

The corresponding diagram is shown in Fig. 2(b). Note 
that that ep L is similar to er with the q' R (l, 3) — ► Uz,(l, 3) 
and u_l(2) — > c^l(2) operations replaced by the (7^(1, 3) — ► 
di(l,3) and <7z,(2) — > gi,(2) operations, respectively. The 
equality of ep R and ep L can be understood by inter- 
changing q <-» gr"jn B'(q' R q L q' R )-B(q' L b)-B(d L q L q' R ) 
and B'(^ 9fl q'[)-B{q' L b)-B{d L q R q'[). The coefficients 
for the | AS 1 1 = 1 case can be obtained by the suitable 



replacement of dr, — > Si, in the B content in Eqs. (A4, 
A5). 

As shown Table IX the octet-octet, octet-decuplet 
and decuplet-decuplet systems are related asymptoti- 
cally [36]. With 

T {,) * l -V ub V: d{s) F T u R v L , 
P (l) « \v tb V; d(s) (F Pl + k^F Pr )u r v l , (A6) 

where as the ratio of the corresponding Wilson coeffi- 
cients of the (db) V - A {qq)v±A {{ds) V - A {qq)v±A) opera- 
tors, we obtain the asymptotic relations, Eqs. (32-37), by 
comparing Table IX with the corresponding amplitudes 
shown in Tables I-V. 



[1] S. Anderson et al. [CLEO Collaboration], Phys. Rev. 

Lett. 86, 2732 (2001). 
[2] K. Abe et al. [Belle Collaboration], Phys. Rev. Lett. 89, 

151802 (2002). 

[3] K. Abe et al. [Belle Collaboration], Phys. Rev. Lett. 88, 
181803 (2002). 

[4] M.Z. Wang [Belle Colaboration] , arXiv:hep-ex/0305060. 
[5] M.Z. Wang et al. [Belle Collaboration], Phys. Rev. Lett. 

90, 201802 (2003). 
[6] N. Gabyshev and H. Kichimi et al. [Belle Collaboration], 

hep-ex/0212052. 
[7] K. Abe et al. [Belle Collaboration], Phys. Rev. D 65, 

091103 (2002). 

[8] K. Hagiwara et al. [Particle Data Group Collaboration], 

Phys. Rev. D 66, 010001 (2002). 
[9] W.S. Hou and A. Soni, Phys. Rev. Lett. 86, 4247 (2001). 
[10] C.K. Chua, W.S. Hou and S.Y. Tsai, Phys. Rev. D 65, 

034003 (2002). 

[11] C.K. Chua, W.S. Hou and S.Y. Tsai, Phys. Lett. B 528, 
233 (2002). 

[12] H.Y. Cheng and K.C. Yang, Phys. Rev. D 66, 014020 
(2002). 

[13] H.Y. Cheng and K.C. Yang, Phys. Rev. D 66, 094009 
(2002). 

[14] C.K. Chua, W.S. Hou and S.Y. Tsai, Phys. Rev. D 66, 
054004 (2002). 

[15] C.K. Chua and W.S. Hou, Eur. Phys. J. C 29, 27 (2003). 
[16] J. L. Rosner, Phys. Rev. D 68, 014004 (2003). 
[17] S. Arunagiri and C.Q. Geng, arXiv:hep-ph/0305268. 
[18] N. G. Deshpande, J. Trampetic and A. Soni, Mod. Phys. 

Lett. 3A, 749 (1988). 
[19] M. Jarfi, O. Lazrak, A. Le Yaouanc, L. Oliver, O. Pene 

and J. C. Raynal, Phys. Rev. D 43, 1599 (1991). 
[20] H.Y. Cheng and K.C. Yang, Phys. Rev. D 65, 054028 

(2002) [Erratum-ibid. D 65, 099901 (2002)]. 
[21] V. L. Chernyak and I. R. Zhitnitsky, Nucl. Phys. B 345, 

137 (1990). 

[22] P. Ball and H. G. Dosch, Z. Phys. C 51, 445 (1991). 
[23] C.H. Chang and W.S. Hou, Eur. Phys. J. C 23, 691 



(2002). 

[24] M. Gronau and J. L. Rosner, Phys. Rev. D 37, 688 

(1988) . 

[25] M. J. Savage and M. B. Wise, Nucl. Phys. B 326, 15 

(1989) . 

[26] X.G. He, B. H. McKellar and D.d. Wu, Phys. Rev. D 41, 
2141 (1990). 

[27] S. M. Sheikholeslami and M. P. Khanna, Phys. Rev. D 
44, 770 (1991). 

[28] S. M. Sheikholeslami, G. K. Sidana and M. P. Khanna, 

Int. J. Mod. Phys. A 7, 1111 (1992). 
[29] Z. Luo and J. L. Rosner, Phys. Rev. D 67, 094017 (2003). 
[30] D. Zeppenfeld, Z. Phys. C 8, 77 (1981). 
[31] L.L. Chau and H.Y. Cheng, Phys. Rev. D 36, 137 (1987). 
[32] L.L. Chau, H.Y. Cheng, W.K. Sze, H. Yao and B. Tseng, 

Phys. Rev. D 43, 2176 (1991) [Erratum-ibid. D 58, 

019902 (1998)]. 

[33] M. Gronau, O. F. Hernandez, D. London and J. L. Ros- 
ner, Phys. Rev. D 50, 4529 (1994). 

[34] M. Gronau, O. F. Hernandez, D. London and J. L. Ros- 
ner, Phys. Rev. D 52, 6374 (1995). 

[35] H.Y. Cheng, Phys. Rev. D 65, 094012 (2002); C.W. Chi- 
ang and J. L. Rosner, Phys. Rev. D 67, 074013 (2003). 

[36] S.J. Brodsky, G.P Lepage and S.A. Zaidi, Phys. Rev. D 
23, 1152 (1981). 

[37] G.P. Lepage and S.J. Brodsky, Phys. Rev. Lett. 43, 545 
(1979) [Erratum-ibid. 43, 1625 (1979)]. 

[38] M. J. Savage and M. B. Wise, Phys. Rev. D 39, 3346 
(1989) [Erratum-ibid. D 40, 3127 (1989)]. 

[39] T.D. Lee, Contemp. Concepts Phys. 1, 1 (1981); 
H. Georgi, Weak Interactions And Modern Particle The- 
ory, Benjamin/Cummings, 1984. 

[40] T. Moroi, arXiv:hep-ph/9503210; S. Deser, V. Pascalutsa 
and A. Waldron, Phys. Rev. D 62, 105031 (2000). 

[41] B. Aubcrt et al. [BABAR Collaboration], arXiv:hep- 
ex/0303028; A. Bornheim et al. [CLEO Collaboration], 
arXiv:hep-ex/0302026; T. Tomura, [Belle Colaboration], 
arXiv:hep-ex/0305036. 



